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Abstract 



An extension of the RINAR(l) process for modelling discrete-time dependent counting 
processes is considered. The model RINAR(p) investigated here is a direct and natural exten- 
sion of the real AR(p) model. Compared to classical INAR(p) models based on the thinning 
operator, the new models have several advantages: simple innovation structure ; autoregressive 
coefficients with arbitrary signs ; possible negative values for time series ; possible negative 
values for the autocorrelation function. The conditions for the stationarity and ergodicity, of 
the RINAR(p) model, are given. For parameter estimation, we consider the least squares 
estimator and we prove its consistency under suitable identifiability condition. Simulation 
experiments as well as analysis of real data sets are carried out to assess the performance of 
the model. 

Keywords: Integer-valued time series, INAR models, rounding operator, RINAR(l) model, 
RINAR(p) model, least squares estimator. 

AMS classification code: Primary 62A/10, secondary 62M20. 

1 Introduction. 

The integer-valued autoregressive, INAR(l), process is introduced by Al-osh & Alzaid 2 (1987) 
to model non-negative integer-valued phenomena that evolve in time. These models are based on 
the binomial thinning operator, denoted o, see [16) . 
The INAR(l) process is defined by 



Here, the so-called counting series ifc) are independent and identically distributed (i.i.d.) Bernoulli 
random variables with success probability oi G [0, 1] and (et) is a sequence of i.i.d. non-negative 
integer-valued random variables and independent of the counting series. Thus, ai o Xt-i is a 
binomial random variable with oi and Xt-i as parameters, ai o Xt-i B{Xt~i, ai). 

The INAR(p) is an analogue of equation ([T]) with p lags. An INAR(p) process is recursively 
defined by 



Xt^aioXt-i+et, VteN, 



(1) 



where 



Xt-i 

ai oXt-i = ^ fifc. 

fc=i 



Xt ^ aio Xt-i -I- 02 o Xt-2 + . .. + apO Xt-p + et, V i £ N, 



(2) 



where, for « = 1, . . . ,p. 




k=l 
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Here (fifc), i £ {I, . . . ,p} and fc G N \ {0} are independent Bernoulli-distributed variables, where 
^ifc has success probability G [0, 1] and (et) is a sequence of i.i.d. non-negative integer- valued 
random variables and independent of all the counting series. 

The general INAR(p) processes where first introduced by Al-osh & Alzaid [3] (1990) but Du & 
Li [7] (1991) proposed a different specification. In the specification of Du & Li (1991), denoted 
INAR-DL, the autocorrelation structure of an INAR(p) process is the same as that of an AR(p) 
process, whereas it corresponds to the one of an ARMA(p,p-l) process in the specification of 
Al-osh & Alzaid (1990), denoted INAR-AA. 
In particular, an INAR(2) process follows the equation 

Xt = ai o Xt-i + a2 o Xt-2 + £t, (3) 

where 

Xt-i 

ai o Xt-i ^J2^ik-~- B{Xt-i,ai) (4) 

k=l 

and 

Xt-2 

02 O Xt-2 = ^ 6fe B{Xt-2, 0,2). (5) 

fe=l 

Note that each variable X^ is thinned twice : ai o A„ for X^+i and 02 o A„ for Xu+2- 

We distinguish two different models according to the construction of these thinning operations: 

The INAR(2)-AA specification introduced by Al-osh & Alzaid [B] (1990). Here, the counting 
series are chosen such that for each A„, the vector (oi o A„, a2 o A„, A„ — oi o A„ — 02 o A„) follows 
a trinomial distribution with parameters (Af_j; ai, 02, 1 — oi — 02). An important consequence 
arises from this choice of the counting series: a moving average structure is induced and the 
autocorrelation function of the process is similar to that of a ARMA(2,1) process. 

Du & Li [7J (1991) propose a modified specification, denoted INAR(2)-DL. The counting 
series are chosen such that {ai o A„, a2 o A„, u G Z} are independent. The correlation properties 
of this process are identical to the AR(2) model. 

Estimators of the parameters of INAR(p) are provided by several authors. For p — I and 
under the assumption that the sequence {et) has a Poisson distribution Franke & Seligmann [Sj 
(1993) analysed maximum likelihood. Du & Li [7] (1991) derived the limit distribution of the OLS 
estimator of a = {ai, . . . ,ap). Brannaas and Hellstrom [S] (2001) considered GMM estimation, 
Silva and Oliveira [TS] (2005) proposed a frequency domain-based estimator of a. Drost & al. 
[S] (2008) provided an efficient estimator of the parameters, and in particular, showed that the 
INAR(p) model has the Local Asymptotic Normality property. 

Thus, the one-step ahead forecast, based on the conditional expectations, for the previous 
models is beset by the problem that forecast values obtained will be real rather than integer-valued 
in all but very rare cases. A mapping into the discrete support of the series is readily obtained by 
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applying Gaussian brackets (integer part of), or by rounding to the nearest integer; the latter will 
be employed along this paper. 

The class of INAR(p) models has several drawbacks. Their innovation structure is complex, 
depending not only on the noise process (et), but also on the counting variables. The autoregressive 
coefficients a^'s are restricted to the interval [0, 1]. In the case of an INAR(l) model for example, 
this restriction excludes negative autocorrelations. 

We introduce the following model 

p 

= ajXt-j +\)+eu te N, (6) 

where (•) represents the rounding operator to the nearest integer, (et) is a sequence of centered 
i.i.d. integer-valued random variables, defined on a probability space (fi. A, P), A and (aj) are real 
parameters. We call this model RINAR(p) (for rounded integer- valued autoregression) . 

RINAR(p) has many advantages compared to the previous INAR models. Its innovation 
structure is simple, generated only by the noise (st). Its one-step ahead least squares predictor is 
given by 

p 

Xt+i = E [Xt+i I Xs, s<T)^ a,Xt_, + A), (7) 

which is integer-valued by construction. We shall also see that the RINAR(p) model can produce 
autocorrelation functions as rich as those of real AR(p), including negatives autocorrelations. 
Moreover, by construction the RINAR(p) model can analyze a time series with negative values, 
a situation not covered by any INAR model. Lastly, it is worth mentioning that the (•) is in fact 
the natural operation often used in the collection of integer-valued series. 

In this paper, we study in details the RINAR(p) model. First in section^ we give conditions 
ensuring the stationarity and the ergodicity of the model. Next in section[3l we introduce the least 
squares estimator for the estimation of the model parameters. This estimator is proved consistent 
under suitable conditions on the model. Because of the discontinuity of the rounding operator, 
particular care is needed for both the formulation of the model identifiability condition and the 
computation of the least squares estimator. A specific algorithm for the last problem is introduced 
in section [H We then present a small simulation experiment in section [S] to assess the performance 
of the estimator. In section[Hl we analyze a well-known time series with RINAR(p) models where 
classical integer- valued models are unsuccessful. Finally, section[7]collects the proofs of all theorical 
results. 

1.1 Some notations 

The following notations and properties will be used along this paper. First, let us define 
N = {x e Z : a; > 0}, Z>o = {a; e Z : x > 0} and = {.T e Z : X < 0}. 

Now, we introduce several useful properties of the rounding operator (•). Note that (a) is clearly 
defined anywhere, unless ii a = k + ^ where A; G Z. By convention, we take {k + ^) — k + l,k 
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and (fc — i) = fc — 1, A; e iT . Note that a (a) is an odd function. 

Let {a} be the fractional part of a G K, {a} G [0; 1[. Here, the fractional part of a negative number 
is a positive as {a} = {—a} = {\a\}, for example { — 1.23} = {1-23} = 0.23. 
Let s be the sign function defined by s(a) = 1 if a > 0, and s(a) = —1 if a < 0. 
Then, for all a G K, we have: 



a = (a) + s(a) {a} ~ s{a)l^^y^i. 

Let [a] be the integer part of a G K, for example [2.8] 2 and [—6.3] = —6. 
Then, for all a G M, we have: 

a = [a] + s{a) {a} . 

The following lemma will be useful and it is a direct consequence of equations 
Lemma 1. Let a; G M and a,b > 0. 

1. \{x) ~x\<\. 

2. |(x)| = (|x|)<|a;| + i. 

3. {a + b)=c+ {{a} + {b}), where c = (a) + (b) - 1|,}> i - l{b}> i - 

4. {a + b} = {{a} + {b}}. 

5. (a)-[a] + (W). 



and 



(8) 



(9) 



2 Ergodicity and stationarity of the RINAR(p) process. 

The study of the RINAR(p) process can be carried out though the following vectorized process 



Yt = 



Xt- 



Xt-i 



(10) 



\ Xt-p+i j \ Xt-p^\ ) 

The process (It) formes an homogeneous Markov chain with state space E = JJ" and transition 
probability function 

p 

7r(a;,y) = P(ei = yi - (^ajXj + A)) ^y^=xu-,v^=x^-i > ^ x ^ {xj),y ^ [yj) E. (11) 

The following proposition gives the conditions which ensure the ergodicity and the stationarity of 
the RINAR(p) process. For x = (xi, . . . ^Xp) G MP , let ||a;||i = + . . . + \xp\. For any measure 
/X and function g on i?, we set /i(g) = J g{x)dfi{x). 



Proposition 1. Suppose that: 



5 



1. The Markov chain (Ft) is irreducible. 

2. For some k>l, E|£t|'= < +00. 

3- E?=iKI<i- 

Then 

1. The RINAR(p) process {Yt) has an unique invariant probability measure fi which has a 
moment of order k (i.e. m(||-||i) < 00^. 

2. For ally & E and f G L^{l-i) we have 

n 
k=l 

where Fy denotes the conditional probability F{. \Yo = y). 



a.s. 



3 Estimation of parameters 

Let 6 = {ai, . . . ,ap, X) € M^+^. In this section, it is assumed that 6 belongs to a compact 
parameters space O, subset of ] — 1, if x R. Let x = {xi, . . . ,Xp) £ W. By convention, we round 
the vector x coordinate-wisely, i.e. (a;) = ((a;i), . . . , {xp)). We note 

p 

f{x;0) = f{xi,...,Xp;e) = (^ajXj + A). 
Then, the RINAR(p) model can be written in the following form 



Y, 



Xt-i 



\ Xt-p+i J \0 J 



where 



M = 



ai . . .ap 
Ip-i 





V J 



and r]t = 





V J 



Let . . . ,Xo, . . . , Xn be observations from the RINAR(p) process. For the estimation 

of the parameter 9, we consider the least squares estimator defined by 



where 



MO) = -Y.{Xt- m.,;0)f = - Y^m - nyt-i; o)h)' 

n ^ — ^ n ^ — ^ 



(12) 



(13) 
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Some notations are necessary. To the norm || • ||i on MP, we associate a norm on (Rp)'^ by setting 
kl •= llyi||i + II2/2II1 for z := {yi, 1/2) in (K^)^, where yi and j/2 are vectors in R^. The actual value 
of the parameter is denoted 9o = (a*, . . . , a*, X*) and Vg^ stands for the probability distribution of 
the chain (Yt) under the actual model. Moreover, any convergence means an a.s. convergence 
under Peo,xi which hold independently of the initial state x. Let us make the following assumptions. 

Assumption [H] 

1. Under Pg,,, the Markov chain (Yt) is irreducible, 

2. for some k>2, E\et\'' < +00, 

3. \ < 1' where a* is the actual value of aj, 

4. the parametric space 6, is a compact subset of ] — 1, x M. 

Assume that [H] holds. Therefore, under Pg^ and from Proposition [Tl (Yt) has an unique invariant 
measure such that Meodl-lli) < 0° for some k > 2. Moreover, the double chain (Zt) with 
Zt — {Yt,Yt-i) has similar properties. So, the chain (Zt) has also an unique invariant measure 
denoted ag^. As Meodl-lli) < follows that o'0q(M*') < 00 for the |.| norm on defined above. 

Let the functions 

9) = (||y - Fix; 9)\\f, z = (x, y) e E\ 9 ^ 6, (14) 

and 

K{e) = ag„g{.;9), 9ee. (15) 

The following proposition give us the limit of the least squares estimating function where the 
convergence holds uniformaly on O. 

Proposition 2. Assume that [H] holds. Then, for any 6 ^ Q, we have 

1. ^^{9) ^ K{9); 

2. K{9) - K{9o) = Meo (/(■; 0) - f{.;9o)f. 
Moreover, 

snp\ipJ9)~K{9)\^0. 

The proofs of this propositions and of all forthcoming results are postponed to section [T] 
Next, we consider the consistency problem of the least squares estimator. As in the RINAR(l) 
case ( see [H]), the identifiability of RINAR(p) model has a non-standard behavior. Because of 
rounding operations, the model identifiability depends whether the autoregressive coefficients a* 
are rational or not. For two parameter vectors 9 — (ai, . . . , a^. A), 9 — {ai, . . . , Op, A ) of O, we 
define their distance to be 

d(0,0') =max{|a, -41,1 < J <p,|A-A'|}. (16) 
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3.1 Strong consistency of the least squares estimator when at least one 
of a* is irrational 

The following proposition adresses precisely the question of identifiability of the parameters of 
RINAR(p) for the case where there exists at least j G {1, . . . ,p} such that a* E R\Q. 
We recall, the following function 
p 

f{x; 9) = ctjXj + A), V a:: = {xi, . . . , Xp) G E and 9 = (ai, . . . , dp, A) S 9. 
i=i 

Proposition 3. Assume that [H] holds. If there exists at least one j € {1, ■ ■ ■ ,p} such that a* is 
irrational. Then, 

!{x-9)^ f{x-9a), •^xeE^9 = 9o. 

Proposition 4. Assume that [H] holds. If there exists at least one j G {1, ■ ■ ■ ,p} such that a* is 
irrational. Then, for all (sufficiently small) e > 0, we have 

inf \K{9)^Ki9o)\>0, 

where Qe ^ {0 : d{9,9o) > e}. 
Theorem 1. Assume that 

1. there exists at least one j G {1, . . . ,p} such that a* is an irrational number; 

2. under Vgg, the Markov chain (Yt) is irreducible; 

3. for some k>2, E\st\'' < +oo; 

4- J2^=i 1*^11 < 1' where a* is the actual value of aj; 

5. the parametric space O, is a compact subset o/] — 1, x R. 

Therefore, the least squares estimators are strongly consistent estimators of the actual values of the 
parameters, i.e. 

9n 6*0, - a.s. 

3.2 Strong consistency of the least squares estimator when all the a* are 
rationals 

First, we recall the main result of the identifiability problem for RINAR(l) model [12], defined 

by 

Xt - +X)+St = f{Xt-i;9) + St. 

Let 6*0 = (ao, Aq) be the actual value of the model. 
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If — ^ , where m G q E Z>o, m and q are taken to be coprime, then 

f{x; 0) — f(x] 6'o), y X G E <=^ a = ao and A 6 Iq. 

Iq is a i-length interval or ^-length interval where Aq S /o- 

The length of Iq depends on the parity of q and the position of {Ao} : 

• If g is even, then Iq is a ^-length interval. 

• If 5 is odd, then we distinguish two sub-cases : 
- If {Ao} e 0; ^ U 1 , then Iq is a ^-length interval. 



2q 



2q 



If {Ao} G , then Iq is a ^-length interval. 



2q ' 2q 

For the RINAR(p) model, we suppose that, for j = 1, ... ,p we have a* = p-, where aj € Z 
and bj G N*, 

aj and 6j are taken to be coprime [aj A bj = 1). Let 0q = (a*, . . . , a*, A*) be the actual value of 
the RINAR(p) model defined by 

Next, we will show that RINAR(p) is equivalent to a RINAR(l) model with a rational parameter 
denoted uq. 

Let y = {xi , . . . , Xp) G E. Thus, 

where 



Ai^fli [] (18) 

From the Bezout theorem, we get 

AiZH + ApZ = dZ, (19) 

where d = Ai A • • • A is the P.G.C.D. of Ai, . . . , Ap. Then, there exists a; G Z such that 

p 

^ a*Xi = VqX, (20) 

i=l 

with 

;.o = Ytp^ (21) 

We note that the numerator and denumerator of vq are not necessary coprime. 
Then, we rewrite vq with its irrcductiblc fraction form 

^0 = ^, (22) 

where a G Z, 6 G N* and a A 6 = 1. 
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Proposition 5. Assume that [H] holds. Let 9 = (a*, . . . , a*, A) G 9. It follows 

yyeE, f{y- 0) = f{y; ^o) ^ V a: £ Z, {i^qx + A) = {vox + A*). 

Proposition 6. Assume that [H] holds. If for all j = 1, . . . ,p we have a* = where aj G Z and 
bj e W, then 

f{x; 0) — f{x; 0o), \/ X £ E aj — a*, \/ j ~ \, . . . ,p and A G /q. 
Jq is a j^-length interval or -^-length interval where A* G /q- 

The length of /q depends on the parity of h and the position of {A*} : 

• If 6 is even, then /q is a ^-length interval. 

• If 6 is odd, then we distinguish two sub-cases : 

- If {A*} G [O; ^[ U [^Ir-; l[' then /o is a ^-length interval. 

— If {A*} G [55; [, then Jq is a ^-length interval. 

Next, in the order to simplify the presentation, we give full details only in the case where b is 
even. 

Let ao = ■ • ■ , Q!p) • From the above discussions, there exists fco G {0, 1, . . . , 5 — 1} such that 
We define, 

Jo = {A : (A) = (A*> and {A} G io} and = {ao} x Jo. (23) 
For example, we take the RINAR(4) model with 

1 3 2^ 4 -i 

"* " 25' ^ 8' " 5 " ^4 ^ 
It is follows, from ([TH| . that 

4 

Jl 6j = 800, = 480, A2 = 1500, A-s = 800 and ^ = -1000. 

So, we have d ~ A\ f\ A2 f\ A-^ t\ A^ — 20. Then, we get, from (|2ip. that t'o = and its irreductible 
fraction form 



Therefore, & = 40 is even, we have 
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Proposition 7. Assume that [H] holds. If ot* = aj G Z and bj e N*, aj and bj are coprime 
y j — I, . . . ,p and b the denumerator of is even, then for all (sufficiently small) e > 0, 
we have 

inf \K{0)-K{9o)\>O, 
eee° 

where 6° = {61 : d{6,Eo) > e}. 
Theorem 2. Assume that 

1. a* — , o,j £ ^ bj G N*, Oj and bj are coprime V j = 1, . . . ,p, 

2. b the denumerator of i2S^) is even. 

3. under Fg^^, the Markov chain (Yt) is irreducible, 

4. for some k>2, E\et\'' < +00, 

5. I'^jl 1' where a* is the actual value of aj, 

6. the parametric space Q, is a compact subset 0/] — 1, 1^ x R. 

Then we get d{9n, Eq) 0, Fg„ — a.s. In other words, ctn is strongly consistent while A„ converges 
to an interval of length i containing A* . 

In the case where b is odd, Theorem [2] still holds where Iq (therefore Eq) is remplaced by the 
corresponding intervals as mentioned in Proposition [SI 

We know that, from [T^], for RINAR(l) model with ao = ^ such that oq and 60 are coprime, 
the length of interval Iq is equal to There is a natural question : 
can we reduce, with p > 1 parameters, the length of the interval Iq ? 
Let us consider the RINAR(2) model with 

a, — — and an = — . 

Then, we get 

Uq = — — = — where d = 0162 A a2&i and a A 6 = 1. 
O1O2 

Let A be an interval in R, we note by \A\ the length of A. We have |/o| — -j. 
Thus 

b/a{bib2) and aA6=l then 6/6162- 

It follows that 

6162 

The RINAR(l) model with the parameter uq = al (resp. 02) produces the interval noted Ai 
(resp. A2). Then, = ^ and 1^21 = ^- Let S = Ai n A2. We have 

/nCB=^|/o|<|B|<min<i-^,l 

61 62 
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Finally, 

_L<l,„l<,„i„{J_.^}, 

Therefore, by increasing the autoregression order p, we can actually expect a reduction of the 
interval Iq where the parameter A remains unidentifiable. 

4 A numerical method to calculate 9n 

To calculate the least squares estimator, we generalize the minimization algorithm proposed in 
[12]. As in RINAR(l) model, the initialization step is given by the Yule- Walker's method. 
Therefore, the generalized algorithm continues through successive dichotomous search steps. 

First, we will explain, in a general way, the algorithm for tpn with a scalar 9. 
The aim is to find 

0n = arg min ipn{0). 
ee[a,fc] 

Here, the initial search interval is [a, b] (i.e. left — a and right — b). The objective function ipn is 
evaluated for every step of the search at three different points: 6k and their middcl left {mid — left) 
and middel right {mid — right) points. According to the minimum value of ipn{Ok), ^n{'mid — left) 
and Lpn{mid — right) ^ 9k, left and right change their actual value (see Figure For example, if 
Vn{ctk) is the minimum then left takes the value of the mid — left and right takes the value of 
the mid — right. This process stops when range = \rigth — left\ < 0.001. 
The following pseudo-code defines the used dichotomous search of 9. 

left <— a; right <— b; range <— abs ( b — a ); 
mid— left <— ( left + previous_theta ) / 2; 
mid— right <— ( right + previous_theta ) / 2; 
while range > 0.001 do 
begin 

V_l <- \Phi_n ( \theta ) ; 

V.2 <- \Phi_n ( mid-left ); 

V.3 <- \Phi_n ( mid-right ); 

i <- j such that V_j is min (V_l , V_2 , V_3) 

case i of 

1 : begin left <— mid— left ; right <— mid— right end; 

2 : begin right <— \theta; \theta <— mid— left end; 

3 : begin left <— \theta; \theta <— mid— right end; 
range <— abs ( right — left ) 

end ; 
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ai a 9i di bi 
initialization [ • • • ] 



first case e • 3 



second case [ » j 



:ai \e, bi 
third case [ • j 



Figure 1: Dichotomous search for a scalar parameter 9. 

4.1 Initialization 

As for RINAR(l), we propose, the initial value 6o = {aifl, . . . ,Q!p,o, Aq) is defined as the Yule- 
Walker estimator as in a real AR(p) model : 

^ "1,0 \ 

p 

; ^Rp^Pp and Aq = X„(l - ^ dj^o), 

where Rp = \p{i — is the sample auto-correlation matrix, pp = (/5(1), . . . , p{p)) and X„ the 

sample mean. 

4.2 Successive dichotomy seeirch steps 

Now, the transition from (di^fc, . . . , ctp^k, Afe) to {ai,k+i, ctp,k+iAk+i), for fc = 0, 1, 2, . . . 
is done in p + 1 phases. The first p phases represent the passage of aj^k to aj^k+i, 3 = ^i - ■ -iP-, 
where the initial search interval is ]— 1, 1[ (i.e. left = —1 and right = 1). The last one represent 
the passage of Xk to Afc+i where the intial search interval is defined to be Aq ± 5|Ao|, which seems 
large enough in most of situations. In every phase, we use the same algorithm described above, 
keeping whenever the results obtained in the phase before. This is the end of [k + 1)*'' iteration. 
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The search stops when the resuhs from two consecutive iterations are very close. More precisely, 
we stop at the /c*'' iteration if : 

d{§k,9k+i) = max||ai^fc - di,fc+i|, . . . , |dp,fc - dp^fc+i|, \Xk - Afe+i|| < 0.001. (24) 

This stopping criterion is satisfied after few iterations. At the end of the iterations, we get the 
vector On ~ (d^, . . . , d^, A„) that minimizes our objective function (^„(q!i, . . . , a^. A). 

5 A simulation study 

For this simulation study, we consider a RINAR(4) model. The error variable, say ei, is generated 
as El — Zi — Z2 where Zi and Z2 are two i.i.d. Poisson random variables. We simulate, n — 500 
observations, of RINAR(4) model with Oq — (^jfj-gj— 3:jf)j values already used in section [X^ 
After 500 independent replications, the average and the standard deviation of the sequence of the 
estimates 6n obtained, are computed and given in the following table. 



^0 


mean 


s.e. 


al = 0.12 


0.1204 


0.0473 


al = 0.375 


0.3687 


0.0439 


al = 0.2 


0.1932 


0.0425 


at = -0.25 


-0.2472 


0.0454 


A* = 2.5 


2.548 


0.2766 



The histograms of these estimates are displayed in Figure [21 Note that inevitably, the a*'s used for 
simulations are rational. Therefore, from Theorem[2 the consistency of the least squares estimator 
9n has to be judged within some inevitable fluctuations of A„. 

Thus, recall that for ao ~ {al,al,al,at) — (^, |, |, —j) we get 1^0 ~ ^ (see Section [X^ . 
We know that 0„ — > Eq = {ao} x /q where Iq is an interval around A* of length 

6 Analysis of the Fiirth data 

In this section, we consider the 505 counts of pedestrians on a city block observed every 5 seconds 
originally published by Fiirth [lOj , see Figure [3l This data set is well known in the branching 
process literature, see Mills and Seneta [Hj (1989). The counts vary from to 7. The mean of the 
series is 1.59 and its variance 1.51. Figure H] shows the sample ACF and PACF of the data. Note 
that in particular, the second partial autocorrelation is significantly negative. 

6.1 The INAR(2) model of Jung and Termayne. 

For the Fiirth data, Jung and Termayne [TT] fitted a INAR(2) model 

Xt = aio Xt-i + a2 o Xt-2 + 6 + £t, 
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where oi, a2 G [0, 1]. Both the two INAR(2) specifications introduced in section [T] are considered. 
With the INAR(2)-AA specification, they have ai = 0.664, 0.2 = -0.119 and b = 0.723. The 
corresponding estimates for the INAR(2)-DL specification are : 0.808, —0.214 and 0.646. None 
of the specifications used yields a satisfactory model, because the fitted value 0,2 is negative which 
is not allowed in an INAR(2) model. 

The authors [11 also fitted a INMA(2) model. The INMA class, as the INAR class, is 
based on the binomial thinnig operator. Its structure is 

= 61 o st-i + 62 o et-2 + b + et, where 61, 62 e [0, 1] . 

The INMA process has been studied by McKenzie [T^ (1988). Al-osh & Alzaid [I] (1988) have 
proposed other specification of the McKenzie's model. For a comparaison of the two specification, 
we refer Brannas and Hall [4J (2001). Jung & Termayne have used the McKenzie specification. 
Then, the Yule- Walker estimates for the parameters are 61 = 1.008, 62 = 0.961 and b — 0.536. A 
main difficulty here is that the estimate for the first thinning parameter bi is outside the admissible 
parameter space [0, 1]. 

6.2 Fit of a RINAR(2) model. 

Here, by using the software R, we will try to fit a real autoregressive model for the Fiirth data. 
Figure m shows that the first and the second sample autocorrelations, are more significant than the 
others and there exists a cutt-off after lag 2 in the partial autocorrelations. First, we consider a 
AR(2) model 

Xt = aiXt^i + a2Xt-2 + b + Et- 

Then, the Yule- Walker estimates for the parameters are 0.808 for ai (s.e. 0.0434), —0.214 for 02 
(s.e. 0.0435) and 0.646 for b (the same estimates for the INAR(2)-DL). Both ai and 02 are 
significant and the AIC value equals 1328.99. 
Now, we consider a AR(3) model 

Xt = aiXt-i + a2Xt-2 + (13X1-3 + b + Et- 

The Yule-Walker estimates for the parameters are 0.828 for ai (s.e. 0.0442), -0.292 for 02 (s.e. 0.0561), 
0.0986 for 03 (s.e. 0.0443) and 0.578 for b. We note, 03 is not very significant. As the associated 
AIC value 1326.05 is very close to the previous one, we will consider the order p = 2 hereafter. 
Therefore, we propose a RINAR(2) model 

Xt = {aiXt-i + a2Xt-2 + A) -I- et- 

In order examine forecast results for this model, we reserved the 400 initial observations as a 
learning set (i.e., to estimate the parameters) and the 105 latest observations as a test set for 
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forecasting. By the algorithm of section^ we find 0„ = (di.n, Q:2.ri, A„) = (0.818, —0.23,0.697). 
The one-step least squares ahead forecast Xt+i of Xt+i equals here : 

Xt+1 = {q-iXt + a2XT-i + A). (25) 

The forecast errors Et+i = Xt+i ~ ^t+i, 400 < T < 504 are displayed on Figure [51 
The mean absolute error equals 

504 

MAE=— J2 |eT+i| = 0.743. (26) 

T=400 

This MAE value is not fully satisfactory indicating that the other RINAR(p) models could 
be considered. Yet, the use of a RINAR(2) model turns to be more natural than the previously 
prposed INAR(2) or INMA(2) fit. In particular, the estimated model is consistent with regard 
to the domain of definition of the parameters which was not the case for the previous INAR(2) 
or INMA(2) models. 

7 Proofs 



Proof of Proposition [T] 

First, we define on MP the functions (p and V with 



vix)^{J2\^j\)''-Ml 



Since V is positive and lim V{x) = oo, V is therefore a Lyapunov fimction ( see Duflo [5]). 

||£c||i— >00 

Let " <" be a partial order relation over R^* defined by 

X <y Xj < yj, Vj = 1, . . . ,p. 



for X = (xj) and y — (yj). We have 



then 



where 



/' -j^ 

J2a,Xt-j+\}+et \ < J2 \^jUt-j\ + |A| + let I + -, 



A 



V {Yt)<A^ iYt-i) + Cu 



|ai| . . . \ap 
Ip-i 



and Ct 



( \M + \et\ + h \ 



\ / 



16 



with Ip-i is the identity matrix of size p — 1. Note that ||CtI|i = \M + k*l + 2' 
Iterating this estimate, we get 

^ (Yt+n) < v4'V (Yt) + A''-\t + ... + ACt + Ct, 
as II • 111 is an increasing function over (Rp, <), 

11^ m+„)l|i < P'V (Yt) + A"-iCt + ... + ACt + Ct\\i 

< ||A> (roili + ll^"-'Olli + • • • + ll^Ctlli + llCtlli 

< ||A> {Yt)h + Sn{A)\\Ct\\i, 

where 

5„(A) =111 III1+...+ III Allli+1, 
with III • I 111 is the operator norm associate to || • ||i. As a result 



ViYt+r.) < ( ||A> (y^lli + SniAm + \et\ + -) 



Then 



k\ k 



{{n-V){x)y < (y,)lli + ^n(A)(|A| 

< (E,||AV (yi)||^)*+^„(A)& 
= IIA'V {x)\\i+Sn{A)b 
<\\\ ^" nil \\x\\i + Sn{A)b, 

where b = (E(|A| + |e4| + i)'')^. We have E|ef 1*= < 00 then < 6 < 00. 

We note that for any matricial norm we have ||| A" |||" — > p{A) as n — > 00, where p{A) is the 
spectral radius of A. As I < 1' then p{A) < 1. It follows, there exists uq such that V n > no 

we have ||| A"" |||= a < 1. So, from Cauchy's criterion we get ^ ||| A" |||i is convergent. 
Therefore, we have 

((7r"F)(a;))*<a||x||i+/3=||x||ifa +. ^ 



Ml 



< a 



So, we get 

n^) - V" ' ll^lli, 

where a < 1 and < /? < 00. Therefore, (Yt) verifies the Lyapunov criterion, with V{x) = \\x\\i as 
Lyapunov function (see Duflo [8 , proposition 2.1.6) and we have that (Yt) is irreducible. It follows 
that (Yt) is positive recurrent with an unique probability invariant measure /i with fiV < 00. The 
conclusion Q follows the classical ergodic theorem for Markov chains. □ 

Some recalls are necessary. £■ = is the state space of the Markov chain (Yj) defined by (fTU]) . 
For all a; = (a;i, . . . , Xp) G E, let ||a;||i — J2^=i Nil- We define | • | the normc on by 

Wz={x,y)eE', 1^1 = ||x||i + bill. 
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As the assumption [H] holds, /ie,, and ctq^ ( the respective unique invariant probabihty of the chain 
(Yt) and the double chain Zt — (Yi^ijYi) under the actual model) have both a moment of order 
k >2. Here we recall the following functions: 

p 

f{x; 9) — cijXj + A), V a; = (xi, . . . , Xp) £ E and 9 — [ai, . . . , Up, A) e Q, 
i=i 

F{x]0) = {f{x;9),xi, . . . ,Xp-i) , 

g{z;9) = {\\y-F{x;9)\\if, ^ z = {x,y) , 

In all of the following Proofs, we denoted a generic constant c whose exact value can change during 
the mathematical development. 

The following Lemma will be useful for the proof of Proposition [2l 
Lemma 2. For all 9 ^ Q, we have 

1. {\\F{-9)\\i?^LHpe„). 

2. \g{-e)\eL\ae,) 

Proof of Lemma [2] 

Note that, for all a G K, \{a)\ < \a\ + i. Thus, by using the compactness of 8, for all x = 
{xi,. . .,Xp) e E, we get 

\\Fix;9)h ^ \fix;9)\ + \x,\ + . . . + \xp.i\ 

< (|ai| + l)|a;i| + . . . + (|ap_i| + l)|a;p_i| + |ap||a;p| + |A| + ^ 

Then, 

{\\Fix;9)\\^r<c{l + \\x\\l). 

So, because k >2 and /ieo(|| • ||i) < oo, the first conclusion follow. 
Now, for all z = (x, y) G i?^, we have 

\\y~Fix;9)h < ||y||i + ||F(x; 

<c(l+||x||i + ||y||i) 
= c{l + \z\), 

it follows 

g{z;9)<cil + \z\^). (27) 
So, because k > 2 and agg{\ ■ |''') < cxo, the second conclusion follow. □ 
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Proof of Proposition [2] 

1. First, we identify the limit of the least squares estimating function. Note that the contrast 
function (^„, defined by p^ . equals: 



if, 



1 " 



(28) 



From Lemma[51 we have g{.] 6) € L^{agg). Therefore, by the ergodic theorem for the double 
Markov chain (Zt), we get 

2. Our aim is to prove that the function K(9) satisfies 

For this, we wiU show that ipn {0) ~ (6*0) ^ fJ-Oa (/(•; - f{-;(^o)f- 
The following definitions and notations will be used in the remainder of the proof. 
Let T = {J-'n)n>o be the natural filtration associated to the RINAR(p) process where 
— c(et,0 < t < n) for n > 1, and J-q is the degenerated (t— algebra. If (M„) is a square 
integrable martingale w.r.t. J^, we denote by ([M]„) its increasing process defined by: 

[M]o = 0, [M]„ = [M]„_i + E(|| [M]„ - [M]„_i ||? | Tn-i) for n > 1. 

Let [M]oo — lini[M]„. On {[M]oo < 00}, we have M„ Moo, where Moo is a finite random 
variable. On {[M]oo = 00}, we have M„/[Af]„ ^ (See Duflo [8,, Theorem 1.3.15, p. 20). 
Recall that, under the actual model 6*0, Yt — F{Yt-i; 9q) + rjt- 
We denote AFt_i = F{Yt-i;eo) - F{Yt-i]6). It follows 

An Bn. 



^n{0) - <Pn(6'o) 



(29) 



with: 



An^Y. W^^t-iWl and Bn = 2^(77* | AFt-i), 
t=i t=i 

where (• | •) is the scalar product associate to || • We have that 
/ f{Yt.,;e)-f{Yt.,;eo) \ 



, so \\AFt-,\\l = if{Yt^,;O)-f{Yt^i;0o)y 



It follows that Mr, 



^ - J2t=i'^tif{Yt-i; 



f{Yt-i;do))- From Lemma [H we have 



/(.;0)P G L^{^0g). Then, by the ergodic theorem, we get 



-^n a.s. 

n 



M«o (/(•;^)-/(-;^o))' 



(30) 
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It is simple to verifies that Mn is a square integrable martinale. Its increasing process [M]„ 
is equal to 

n 

[M]„-^(/m_i;0)-/(yt-i;0o))'lE|£tr 
t=i 

As almost surely, 

-[M]„ ((/(.; 0) - /(.;0o))'r) >0, where T^E\et\^<oo, 

n 



it follows that -M„ — > 0, therefore 

Bn 



3. Our aim is to prove that 



sup\ipn{e)^K{e)\ ^0. 



Some notations are necessary. For any z = {x,y) e E"^ , let \z\ ~ \\x\\ + 
Let P„ be empirical measure generated by the observations Zi , . . . , Z„ 



(31) 



1 

n 

1=1 

It follows, from equation (|28p . the contrast function p„ equals: 

Pn(^?)-P„g(.;0). (32) 

From equation (|27p . we get 

|g(z; 6*)! < yl(z), V z = (a;, y) £ andfl e 6, (33) 

where 

A{z) = c(l + |zp), where c is a constant. (34) 
Note that, A G L^{(T0„). Let g > be fixed. It arises, 

|^„(0)-if(0)| = |(P„-a9„)5l 

= |(Pn - cr0o)5l|z|<q + (Pn " (^eo)9^\z\>q\ 

< |(P„ -tTe„).ga|,|<,|+P„(|.g|l|,|>,) + aeo(|.9|l|,|>,) 

< |(P„ -f7eo)5ak|<9l+]Pn(A(z)]li,|>,)+ae„(A(z)]li,|>,). 

We denote = '^eo{z) and = P„(2;). Moreover, 

|(P„-aej5a|^|<,l - I ^ g(z;0)(p?-p.)| 

< ^A(z)b:-p,|. 

This means that 

sup l^niO) - K {9)1 = snpliF^ - ag„)g\ < ^ A(z)|p^ - + (P„ + <7eJ(A(z)l|,|>,). 

|.|<, 
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When n ^ oo, we have — > Pz almost surely, then the finite sum 
J2\z\<qMz)\Pz - Pz\ 0, a.s. Therefore, 

limsupsup |(P„ - CTeo)5| < + 2a0„{A{z)liz\>q), a.s. 
By making q y oo, we get almost surely, 

limsup sup [(P„ — CTSo).?! =0. □ 

n^oo eee 

Recall that O is compact. Thus, for all 9 — {ai, . . . ,ap, X) G 8, there exists two positive 
constants A and B such that |A| < B and \aj\ < A, V j = 1, . . . ,p. {■} represents the fractional 
part operator. For all a G K, {a} G [0; 1[ and {a} = {\a\}. 

Proof of Proposition S] 

Our aim is to prove that if there exists at least one j G {1,2, ... ,p} such that a* G M\Q, then for 
all 6 — (ai, . . . , Up, X) G Q, we have 

p p 

aiXi + A) = a*Xi + A*), V x — (xi) G E =^ \ ~ \* and ai — a*, Vi — 1, . . . ,p. 

i=l 1=1 

The idea is to prove that ii 9 ^ 6q then there exists xq E E such that f{xo; 9) ^ f{xo; 9o). 

We assume, without loss of generality, ai ^ al. Let x = (a;i,0, ... ,0) G E. Thus, by using the 

compactness of 0, we get 

\f{x;9)-f{x;9o)\ > \ai - al\\xi\ - 2B - 1. 

Therefore, if (a;) ^ (a*), then there exists xq £ E such that \f{xo;9) — /(xo;0o)| > 0, this is a 

contradiction. Moreover, if xq = Oe, then f{xQ;9) — /(a^oS^o) implies (A) = (A*). 

In order to simplify the notations, we assume that a\ G K\Q. Let y — {x,Q, ...,{))€ E. 

We assume, without loss of generality, > 0, A > A* > and {A} , {A*} G [O; \ [. 

It follows, A - A* = {A} - {A*} and 

f{y; 9) = /(y; 9^) ^ {{a\x} + {A}) = {{a\x} + {A*}). (35) 

As a\ G it arises {{a\x})xe7. is dense in [0; 1[. If {A} ^ {A*}, then there exists yo E Z such 

that 

{alyo} + {^} > ^ and {a*yo} + {A*} < ^ 
Finally, there exists xq = {yo, 0, . . . ,0) E E such that f{xo; 9) ^ f{xo;9o). □ 
Next, we consider d the distance on the parametric space Q, defined by 

d{9,9') = max||aj - a]\,l <j <p,\X~ A'|| , V 61,6'' G 9. 

Recall that, from Proposition [2 we have 

Ki9) - Ki9o) = fieo (/(•; 0) - /(.; ^o))' ■ (36) 
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Proof of Proposition SI 

Our aim is to prove that if there exists at least one j G {1,2, ... ,p} such that a* E R\Q, then for 
all (sufficiently small) e > 0, we have 

inf \K{e)^K{eo)\>0, 

where = {9 : d{9, 9o) > e}. We need distinguish three situations for the event Qe- 
Thus, we note that 9^ Ti U r2 U Ta, where 

Ti^{e:3ie{l,...,p},\a,-a*\>e}, (37) 

r2 = {9:Vte{l,... ,p} , \a, -a*\< e, \\ - \*\ > e, (A) ^ (A*)} , (38) 

and 

Ts = : e {f, . . . ,p} , \a, - a*\ < e, |A - A*| > e, (A) = (A*)} . (39) 
We are going to prove 

inf \K{9)-K(dQ)\ >0,i = f,2,3. 
The idea of the proof is based on equation . Therefore, the aim is to find xq £ E such that 

\f{xo;9) - /(a;o;6'o)| > 0, uniformly on T^,i = 1,2,3. 

1. We consider the first case, 9 £ Ti. As is compact implies that |A| < B, where B is a 
constant. Let x — {0, . . . ,0,Xi,Q, . . . ,0) £ E, we get 

\fix;9)- f{x;9o)\>e\x,\-2B-l. 

Therefore, there exists yo > such that V |a;i| > i/q, we have \ f{x;9) — f{x;9o)\ > 
uniformly on Fi. It follows, since the support of p-e^ is not bounded, 

inf \K{9)-K{0o)\>O 

Pti 1 

2. For 9 £ Fa, let Oe = {0, . . . ,0) £ E, we have 

|/(Ob;0)-/(Ob;(?o)|-|(A)-(A*)|>1. 
Because, if two real numbers do not have the same integer So that 

inf \K{9)-K{9a)\ > 0. 

3. Finally, we consider the last case, 9 £ T^. Here, we assume without loss of generality, 
al £ K\Q. Let E =^ {x ^ (x,) e £: : xi e Z and = 0, V j = 2, . . . ,p}. 

On E' , the actual RINAR(p) becomes a RINAR(l) model with al £ R\Q. 

Then, we retake the same proof of Proposition 4 from |12] for the similar case, remplacing in 

each step xq by a vector yo — (xq, 0, . . . , 0). □ 
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Proof of Theorem [T] 

The conclusion On Oo almost surely results from Propositions [4] and [2] by standard arguments of 
the theory of M-estimators (see e.g. Van Der Vaart [TZ , Theorem 5.7). □ 

Proof of Proposition [5] 

Let y = {xi, . . . , Xp) £ E. Then, from (fT7|) (fT5|) and we have that for every y £ E there exists 
X € Zi such that X]r=i (^i^i — '^ox- Let x £ 1^. From ()2ip . we have that i/qx — ;ppp — T~i'^ ^) ''^^'-^ 
again from (jl9p we get that there exists y — {xi, . . . , Xp) G E such that d x = AiXi + . . . + ApXp. 
It is follows that 

Vox = YW—-{AiXi + . . . + ApXp). 
Finally, from (fT7|) . we get that for every x £TL there s ?/ = (xi, . . . , Xp) G E such that 

Proof of Proposition [6] 

Note that K{9) ~ K{Oq) means, from Proposition /(.; 0) = /(.; 6*0) /iso ~ ^-^ 

i.e. y y — {xi , . . . ,Xj, . . . , Xp) € E, we have 

p p 

aiXi + A) = a*Xi + A*). 

i=l i=l 

Let = (0, . . . , 0, Xi, 0, . . . , 0) G for i — 1, . . . ,p. Letting \xi\ — > 00 for each we find 

Qfi = a* V i = 1, . . . ,p. 

So, from Proposition [51 we get 

VyeS, /(y;0) = /(y;0o)^VxeZ, {u^x + \) = {u^x + \*) . 

It is follows, from Proposition 5 of [12 , that 

V x e Z, {uox + A) = (z/oa; + A*) ■^=^ A e /o. □ 

Let = (xj, . . . , x^) G E. By Bczout theorem, we know that there exits xo G Z such that 
SiLi "^i^? ~ i^qXq. The following lemma will be useful for the proof of Proposition [71 

Lemma 3. The function a — {ai, . . . ,ap) — > defined on MP, is continuous at 

ao = (a* , . . . , a* ) if Xq ^ mb, with m G N and b is denumerator of irreductible fraction form of vq . 
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Proof of Lemma [3] 

Let x*' — (x?, . . . , Xp) e E fixed. We know that the function h : a — > {a ■ x"} = |X]j=i Q^j^?!! 
defined on R^* and witli values in [0, 1[, is discontinuous at ao, if X]j=i '^j-'^^ ^ ^- From Bczout 
theorem, there exists xq Z fixed, such that 

p 

E* ° 

So, if X{) ^ m6, we get vqXq Z, it foUows h is continuous at ao- D 
Recall that, d the distance on the parametric space 0, is defined by 



d{d,d ) = max||aj - < j < p, |A - A'|| , ^ 9,9 e 



e. 



Proof of Proposition [7] 

We have, — j; where a G Z and 6 G Z>o (fo is even), a and 6 are comprime 
Thus, there exists fcg G {0, 1, . . . , g — 1} such that 

fco fco + 1 



6 ' 6 



{A*} G = 
We recall, 

/o = {A : (A) = (A*), {A} G io} and Eo ^ {9 e Q : oij ^ a* , V j = 1, . . . ,p and A G /o} . 
Our aim is to prove that for all (sufficiently small) e > 0, we have 

inf \K{9)-K{9q)\>0, 

where — {9 : d{9,Eo) > e}. We need distinguish three situations for the event 8^. We note 
that = Ti U U Ta, where 

Ti = {0:3iG {!,..., ri,|a, -«*!>£}, (40) 

= {0 : G {1, . . . ,p} , \a, - a*\ < e, d(A, /o) > e, (A) ^ (A*)} , (41) 

and 

= {^^ : Vz G {1, . . . ,p} , |a. - a*! < e, d(A, /q) > £, (A) = (A*)} . (42) 
We are going to prove 

inf \K(9)-K(9q)\ > O.i = 1,2,3. 
From Proposition we know that 

Ki9) - Ki9o) ^ fieo (/(•; ^) - /(•; ^o))' ■ 
So, the aim is to find xq £ E such that 

\fixo;9)- fixo;9,)\>0, V 0, G T,, i = 1, 2, 3. 
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1. For G Fi U then for the same arguments of the equivalents subsets of the proof of 
Proposition m we have that 

2. For 9 e F3, without loss of generality we assume that 1^0 > 0, A* > 0. Note that 



d{e,Eo)=d{X,Io)= inf IA-A*|. 
We distinguish four cases depending on the position of the fractional part of A* and A : 

(a) Case A : in C [O, ^ [ and A > /q. 

(b) Case B : iq C [O, | [ and A < Iq. 

(c) Case C : io C [5, 1 [ and A > Iq. 

(d) Case D : io C [i, l[ and A < Iq- 

• Case A : if) C [O, ^ [ ; A > /q. For this case, it is simple to verified that 



b > 4, {A} e 



Therefore, 



d{0,Eo)= inf {A}-{A4>e. 



(43) 



There exists n e {ko + 1, . . . , | - l} such that {A} G ['f, ^^ii Let = | - n fixed. 
We know that, from Bezout theorem, for xq used here there exits = {xf, . . . ,Xp) EE 
such that 



E* 



For all {A*} G 2o we have 



E">o| + {A4 = |ixo| + {A4 < i 



(44) 



and from we get 



f:a*x°|+{A}=|i.xo} 



+ {A}> 2+£- 



(45) 



Also, for Xq = ^ — n, from Lemma [U the function a = (ai, . . . , Up) — *■ '^i^j | 

is continuous in ao = (aj, . . . , a*). There exists 77 = ry(£, t'o, A*) < e such that V j = 
1 , . . . , p I — a* I < 77 it arise 



(46) 
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From ([IS]) and (HH]) we get 



E".-" +w = ( E".-? - E"M 



= 1 



Then 



p 

E "J^i " 




p 

-(E 








p 






E "^J^i 






.j=i 






1 P 






iE"^-^°r 


-{A})- 









+ (A) - (A.) 



= 1 



• The other cases : 

— Case B : zo C [O, ^ [ and A < /q. 

— Case C : io C l[ and A > /q. 

— Case D : zo C 1 [ and A < /q. 

can be treated in a similar way, with a suitable choice of xq, and a = (ai, . . . , Up) suffi- 
ciently close to ao = (a^^, . . . , a*). Then, we deduced that there exists 77 = 77(6, i^Oi A*) < 
£ and a constant cq > 0, uniformly on Fg = |6' : 6' G F3, V j G {1, . . . ,p} |aj — a*| < 77}, 
we have \K{6) - K{ea)\ > eo > 0. 

For 6* e F3 = F3 \ F3 : 6* e F3, V j e {1, . . . ,p} |a-, - a*] > 77} the same argument 

for Fi can be applied here and then there exists c?o > such that 



inf \K{e)-K{eQ)\ >do>0. □ 
eer" 



Proof of Theorem [2] 



The conclusion 0„ Eq, almost surely, results from Propositions [2] and [7] by the same arguments 
used in the proof of Theorem [TJ □ 
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Figure 2: Histograms of the parameters estimates aj,j = 1,...,4 and A (left to right, top to 
bottom). 500 repUcations of time series of length 500. Actual values of the parameters are 6o = 

(3_ 3 1 -1 5\ 
V25'8'5' 4 ' 2)- 
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Figure 3: Plot of the Fiirth data. 
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Figure 4: Sample autocorrelation and partial autocorrrelation functions of Fiirth data. 
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The forecast errors of the 1 05 last observations. 
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Figure 5: The forecast errors ix+i = Xt+i — Xt+i, 400 < T < 504, of the Fiirth data, from a 
RINAR(2) fit. 
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